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We propose a method to expand correlation functions with respect to the spatial components 
of external momenta. From the coefficients of the expansion it is possible to extract Lorentz- 
invariant form factors at zero spatial momentum transfer avoiding model dependent extrapolations. 
These objects can be profitably calculated on the lattice. We have explicitly checked the validity of 
the proposed procedure by considering two-point correlators with insertions of the axial current, the 
form factors of the semileptonic decay of pseudoscalar mesons, and the hadronic vacuum polarization 
tensor entering, for example, the lattice calculation of the anomalous magnetic moment of the muon. 



I. INTRODUCTION 

Physical observables are calculated on the lattice by 
extracting hadronic matrix elements from the large eu- 
clidean time behavior of correlation functions. The ma- 
trix elements are Lorentz-covariant quantities, i.e. ten- 
sors of a certain rank and with definite properties under 
discrete symmetries. These are parametrized in terms 
of form factors, the Lorentz-invariant coefficients associ- 
ated with the decomposition on a certain basis. The ba- 
sis is built by considering products of external momenta, 
polarization vectors, metric tensor, etc. The form fac- 
tors multiplying powers of the the external momenta are 
difficult to extract when these variables vanish. 

The class of problems we address in this letter can be 
exemplified as follows 

C»-(p,---)=p»F-(p 2 ,...), 



(1) 



where C M **' is the quantity that can be computed directly 
on the lattice (a correlator or a combination of correla- 
tors), F " is the quantity one is interested in, while the 
dots represent additional arguments that may or may not 
be present and that characterize the given observable. 

Typically it is possible to obtain good numerical preci- 
sion on the form factors for values of p M in a limited range. 
Lattice correlators tend to be extremely noisy for large 
values of the modulus of p^. On the other side, although 
the use of flavour twisted boundary conditions [TJ |2] al- 
lows to overcome the quantization of spatial momenta on 
finite volumes and to compute F " for very small values 
ofp M , the numerical signal for C^'" goes to zero for nearly 
vanishing values of p M and the results in this region are 
again very rough. In particular there are situations in 
which one is interested in the value of F " at p^ = 
that, by using standard lattice techniques, can only be 
obtained by recurring to (model dependent) extrapola- 
tions. 

A concrete example of this kind is offered by the cal- 
culation at zero recoil of all the independent form factors 
parameterizing the semileptonic decays B — »• D^tv [3]- 



Another phenomenologically relevant example is given 
by the hadronic vacuum polarization contribution to the 
muon g — 2 that is extracted from the integrated two- 
point correlator of the electromagnetic currents of the 
quarks, 
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(2) 



In this case the accurate knowledge of 11(0) is extremely 
important because divergent contributions are canceled 
by extracting the physical information from the differ- 
ence II(p 2 ) — 11(0). By using standard approaches, 11(0) 
can only be obtained by extrapolating the results at non 
vanishing p M and by relying on fitting formulae that un- 
avoidably introduce systematic errors (see ref. [8 for a 
recent discussion of this point). 

In this work we discuss a method to be used in or- 
der to solve the class of problems discussed above. The 
method allows to calculate directly on the lattice the nth- 
derivative of a lattice correlator with respect to external 
spatial momenta. By using the notation of the previous 
formulae, we shall describe in the following sections a 
procedure to calculate on the lattice the left hand side of 
(no sum over k) 



dp k 



(3) 



directly at p k = without the need of any extrapolation 
and related ansatz for the fitting function. By applying 
this method it is possible to compute on the lattice the 
vacuum polarization scalar form factor at zero momen- 
tum according to 



n(o) 



d 2 c 12 (p) 



dp 1 dp 2 



(4) 



p 2 =0 
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II. EXPANSION OF THE CORRELATORS 

In order to expand lattice observables in powers of 
external spatial momenta, let us start by considering a 
generic correlator depending upon 

C ^ = ZsXX^ / dU p [U]C[x;U] 

X 

dUP[U]C[p;U] . (5) 



In previous expressions P[U] represents the gauge "prob- 
ability density" and incorporates the fermionic determi- 
nants together with the normalization factor of the func- 
tional integral. Note that P[U] does not depend upon 
external momenta. Our method consists in expanding 
C[p;C7], i.e. the correlator at fixed QCD gauge back- 
ground, with respect to p according to 



C\p;U]=C {0) [U]+p k cf ) [U] 



PhPk r ,(2) 



(6) 



The coefficients a (n) [C7] can be interpreted as the deriva- 
tives of the original correlator with respect to the external 
momenta. Note that this is true also when the external 
momenta are quantized as a consequence of the finite 
volume and of the boundary conditions satisfied by the 
quark and gauge fields. 

After fermionic integration and Wick contractions the 
lattice correlator C[p;U] can be expressed as the product 
of quark propagators and operator insertions, i.e. 



C(p) 



J dUP[U}^^2 e^-^Tr {S[x, y;U]T- 
<^£Tr{S[*,y;[/,AP]r...>. 



(7) 



Her e T represent a generic operator insertion (see sec- 
tion [lV| below) while S[x,y]U] is one of the fermion prop- 
agators entering the expression of the given correlator. 
The fermion propagator S[x,y;U] (S[U] in short) satis- 
fies the following equation 



(8) 



where D[x, y\U\ is the lattice Dirac operator used to dis- 
cretize the quarks action. The propagator S[x, y; U, X p ] 
introduced above (S[U, X p ] in short) is simply given by 



S[x,y;U,X p ]=e ip ^- g) S[x,y;U\. 



(9) 



Dirac equation with the link variables rescaled by a phase 
factor 

U k {x) .— ► \ p k U k (x) = e ip " U k (x) , 

D[x, y; U, \ p ] S[y, z; U, \ p ] = S x , z . (10) 

y 

Note that S[U, X p ] is defined analogously to a lat- 
tice propagator associated with a quark field satisfying 
flavour twisted boundary conditions [TJ [2] , with the dif- 
ference that now p is not a continuous variable. In this 
work we shall use periodic boundary conditions for all 
the quark fields and, consequently, we shall always work 
within a unitary setup and with the spatial momenta 
quantized according to 



P=Y V n l' n 2,n 3 ) . 



(11) 



Expressions for the coefficients C^[U] can now be 
readily obtained by expanding the lattice propagators 
S[U, X p ] and possibly non ultra-local operators entering 
into the original expression of C[p; U] (see the following 
sections for examples) in powers of the external momenta. 



III. EXPANSION OF THE PROPAGATORS 

The expansion of the quark propagators can be ob- 
tained by starting from the expansion of the lattice Dirac 
operator, 



X p k = 1 + ip k ~ 7>Pk 



D[U 1 X p ]=D[U]+p k 



dD 

dpk 



p=0 



2 dpi 



p=0 



(12) 



and by using the identity D[U, \ P ]S[U, \ p ] = 1 leading to 



dS __ 
dp k 

i d 2 s 

2 dpi 



dD 

op k 
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(13) 



In the previous expressions there is no sum over k and 
we have used the following compact notation 



dO = dO["- ;U,X p ] 
dp k ~ dp k 



(14) 



p=0 



The previous definition has been made because it is use- 
ful to obtain S[U, X p ] directly, i.e. by solving the lattice 



For our numerical calculations we have used clover Wil- 
son fermions, corresponding to the specific Dirac opera- 
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tor D[U,X P ] given below 
^(x)D[x,y;U,\ p Mv) 



<K*)(4 + m)^[x) + ^ Yl ^)a^F^(x) ${x) 



1-7^ 



^^ + M )(A^)t[/t (a;) 



1 + 7 A 



(15) 



where Aq = 1, A^ = e* Pfc , F^ v (x) is the clover defini- 
tion of the lattice gluon field strength tensor and csw 
is the Sheikholeslami-Wohlert coefficient [9]. The nth- 
derivatives of this operator are proportional either to the 
"point-split" vector current V k (x,0) (for odd n), or to 
the "tadpole" current T k (x,0) (for even n), 



T dD 

OPk 



dpi 



tp%- — ip 



(x) = iV k (x,6) 
(x) = T k (x,6) 

(x) = -{-if 



V k (x,0) odd n 
T k (x, 0) even n 



(16) 



where the expressions of the two currents are given by 

v»(x,p) = [$rtw>](x,d 

= +^ x + ri(\P^ul{x) l -±J^iP{x) 
-j(x)\^(x)^-^iP(x + p), 

T»(x,p) = [^V] (x,p) 

= +^ (a; + M )(Apt C /t (x) l+^^ (x) 

+4 > (x)\lU ll (x)^^4>(x + (i) . 



(17) 



Graphically, the expansion of the propagator of eqs. ([13]) 
can be represented as the sum of iterated insertions at 
zero momentum of the operators V k (x,0) and T k (x,0) 



according to (no sum over k) 

Pk 



Wk 




(18) 



In the previous expressions we have introduced the fol- 
lowing graphical notation 



y ■ x = S[x,y,U] = (<ip(x)<ip(y)) , 

# = £r*(*,6), 

z 

M = ^r^(z,o), (19) 

z 

that we shall repeatedly use also in the following sections. 
IV. EXPANSION OF THE VERTICES 



As discussed above, the expression of a correlator at 
fixed QCD gauge background involves products of traces 
of propagators S[U, X p ] and operators T, as rewritten here 
below in a very compact notation, 

C(p) = (C\p; U]) = (Tr {S[U, X p ] T p ■ ■ ■ } . 

If the operators are not ultra-local and connect different 
lattice points, they may also carry a dependence upon 
Pk- For this reason we have used in previous expression 
the notation T p . In these situations the calculation of the 
coefficients of the momentum expansion of the correlators 
requires the knowledge of the derivatives with respect to 
external momenta of all p-dependent objects, including 
the vertices, as shown below for the first derivative 



C^[U] =Tr 



,0T 

dpk 



(20) 



One relevant example of p-dependent operator is just 
given by the point-split vector current discussed above. 
It undergoes the following expansion 

V h (x,p) = cos(p k )V k (x,6) -ism(p k )T k (x,6) 



(l-^- + ...)V k (x,6)-i(p k + ...)T k (x,0) 



(21) 



to which, in the case of the integrated insertions of these 
operators, we associate the following graphical represen- 
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tation 



relators 



^2^v(z,P) = # = cos(p fe ) # - isin(p fe ) 



-(ip*)B-f •+••• • (22) 



The previous expressions will be used below in order to 
expand the integrated correlator proportional to the vac- 
uum polarization tensor. 



APPLICATIONS 
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,ip(y-x) r ; 5 



[V»i7 VsRs) [V>27 Vl](tf)> 



L 6 



Tr {%, x; C/] 75 S[x, y; [/, A p ] 75 }> , 



_ ZaZp 
L 6 



<£■ 



3 ip(y-x) 



In this section we check the validity of the method by 
calculating on the lattice the derivatives with respect to 
external spatial momenta of a set of correlation functions. 

We start by considering two-point correlators of a pseu- 
doscalar density and of an axial current. This example 
has no relevance from the phenomenological point of view 
but it is very useful to explain how the method works in 
practice and to check its effectiveness in a simple and con- 
trollable situation. As a second application we consider 
three-point correlation functions and use our method to 
extract at zero recoil both form factors parametrizing 
pseudoscalar-to-pseudoscalar semileptonic transitions. 
Finally we compute the second derivative of the vacuum 
polarization tensor and extract the scalar form factor 
needed for the muon g — 2 calculation directly at zero 
momentum. 

The numerical results have been obtained (on a single 
volume, at fixed lattice spacing and fixed sea quark mass 
in the rif = 2 theory) in order to check the validity of 
the proposed procedure. For this reason, also in the case 
of the phenomenologically relevant applications discussed 
below, we do not quote numbers in physical units neither 
attempt to estimate possible systematic errors. 

For our numerical investigations we have used the 
gauge configurations corresponding to the D2 ensem- 
ble of the CLS initiative. These in turn correspond 
to the simulation parameters V = 48 x 24 3 , /3 = 5.3, 
k sea = 0.13590 and c SW = 1.90952. All the other infor- 
mation concerning this gauge ensemble can be found in 
refs. 501 E] . 



A. Two— point correlation functions 



A simple and illustrative check of our procedure have 
been obtained by considering the following two-point cor- 



ZaZp 
L 6 



<]T Tr {%, x- U] l5 S[x, y; U, A^ 75 }>. 

x,y 

(23) 



For large euclidean time separations t = xq — yo we can 
neglect sub-leading exponentials contributing to the cor- 
relators and write 



C PP {t,p) = 



CP 
2E 



-Et 



FG 
2E 



-Et 



(24) 



where is the momentum of the pseudoscalar meson, 
M its mass, po = E = + p 2 its energy while F and 

G are Lorentz-invariant form factors depending upon the 
meson mass square p 2 = M 2 . 

In this simple situation our method is not needed be- 
cause the decay constant F of the meson can easily be 
extracted from the correlators Cpp(t,0) and C\ P (t, 0) 
both evaluated at vanishing spatial momentum. On the 
other hand, we can check our proposal by computing the 
first coefficient of the momentum expansion of C\ P (t,p), 
namely (no sum over k) 



d 
dpk 



C k AP (t,p) 



p=0 



• Tr | S[y, x ; U] l5 A S[x, y; U] 7 k l5 } ) 



x,y 



dpk 




(25) 



where, according to the formulae derived in the previous 
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FIG. 1: The black points correspond to the effective mass 
obtained from the correlation function Cpp(t, 0) while the red 
points correspond to the mass extracted from the ratio B,Ap(t) 
(see definition in the text) computed at p = 0. The two 
definitions of effective meson mass agree for large euclidean 
time within statistical errors. The error bars are of the same 
size of the symbols used for the plots. Data are in lattice 
units. 



sections, we have 



— S[x,y,U] = —S[x, y] U,\ p ] 
dpk dp k 



p=0 



= -J2s[x,z;U}iT k v (z,6) S[z,y;U] . 

z 

(26) 

We remark that the quantities above have been deter- 
mined in the periodic theory with X p = 1 or, equiva- 
lent ly, p = 0. The large time behavior of dC\ P /dpk is 
now expected to be 
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dpk 



C AP (t,p) 



p=0 



FG 
2M 



-Mt 



(27) 



In order to check this behavior, in FIG. [T] we compare 
the ratio 



R AP (t) 



c° AP (t,o) 



^C k AP (t,p) 



M 



(28) 



p=0 



with the effective mass mpp(t) extracted from the 
pseudoscalar-pseudoscalar correlator Cpp(t, 0) at zero 
spatial momentum. As it can be seen, for large time 
separations there is agreement within the statistical er- 
rors between the two definitions of effective mass, thus 
confirming the theoretical expectations and the validity 
of our method. 

An important remark is in order here. Once a given 
correlator it has been made finite by multiplying it with 
the appropriate renormalization factors, the derivatives 



of the renormalized correlation functions with respect to 
external momenta are also finite, i.e. there is no need of 
specific renormalization prescriptions for the coefficients 
of the momentum expansion. In particular, the ratio 
RAp(t) is finite (renormalization factors identically can- 
cel between the numerator and the denominator) and has 
indeed been plotted in FIG. [I] without multiplying it with 
any factor. 



B. Three— point correlation functions 



In this section we apply our method to three-point 
correlation functions. In particular we consider the form 
factors entering the rate of the semileptonic decay of 
a pseudoscalar meson into another pseudoscalar meson, 
Pj — ^ P F £v. The process is mediated by the vector part 
of the weak V — A current and the corresponding matrix 
element can be parametrized in terms of two independent 
Lorentz- invariant form factors, 

<Pf|V|P/> = (Pi +p F rf+(q 2 ) + (pi- PF rf-(Q 2 ) , 

(29) 

where pj and p p are respectively the initial and final 
momenta of the hadrons and 



(f =Pi ~Pf • 



(30) 



Here and in the following the dependence of form 
factors upon Lorentz invariants is understood, f± = 

f±(p 2 pP F , Q 2 ) with P! = M2 I and P% = M F' 

The form factors are usually extracted by consider- 
ing ratios of three-point correlators evaluated at different 
initial and final masses and momenta. By introducing a 
compact notation for the initial, final and spectator quark 
propagators (rrtx are the quark masses), 

5(0) ^ S[y,x;mi; U] , 

5 J (pj) i y Siz,y;m i; U,\ P1 } , 

S F (p F ) ^ S[x,z;m F ,U,\ PF ] , 
these correlators can be written as follows 
C IF (t,p!,p F ) 



(31) 



_ J_ \^ ip F {z-x)+ip I {y-z) I 



x,y,z 



(32) 
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FIG. 2: The black points correspond to the effective mass ob- 
tained from the correlation function Cpp(t, 0), the red points 
correspond to the mass extracted from the ratio Rap (£) while 
the green points correspond to the mass extracted from the 
ratio Rv(t) (see definition in the text) computed at p = 0. 
The error bars are of the same size of the symbols used to 
plot the data. Data are in lattice units. 



In the previous expression the distance between the in- 
terpolating operators of the initial and final pseudoscalar 
mesons has been set to half of the time extent of the lat- 
tice, i.e. xo — yo = T/2, we have called t = zo — yo, the 
improvement of the vector current has been neglected 
and the multiplicative renormalization factors of the in- 
terpolating pseudoscalar densities and of the local vector 
current have been omitted because these cancel in the 
ratios from which the form factors are extracted (see be- 
low). 

Our gauge ensemble has been generated with two dy- 
namical quarks having the same mass while in eq. (32) 
we are referring to three different quark flavours. 



The 

third quark, in the case of different initial and final meson 
masses, has been treated in the quenched approximation. 

When the masses of the initial and final mesons are 
equal, the conservation of the vector current implies the 
presence of a single form factor 



(Pi\V»\? I ) = ( Pl + PF T U(q 2 ) 



(33) 



exactly normalized to 1 at q 2 = 0. The simple situation of 
equal initial and final meson masses at zero recoil has no 
phenomenological relevance but gives us the opportunity 
of repeating the check already performed in the case of 
two-point correlators also in the case of three-point cor- 
relation functions. To this end, we set the initial spatial 
momentum to zero, p} = 0, and consider the behavior 
of the correlators for large relative time separations by 
neglecting sub-leading exponentials, 



C°j(t,0,p) = (M + E)f + (q 2 ) 



G 2 



AME 



-E(T/2-t)-Mt 



CUtAp) = P k f+(q 2 ) 



G 2 



AME 



-E(T/2-t)-Mt 



(34) 



where M = Mi = Mp, p is the spatial momentum of 
the final meson, E — \J M 2 + p 2 its energy and G the 
matrix element of the interpolating operator. We have 
computed the first coefficient of the momentum expan- 
sion of Cjj(t,0,p), i.e. (no sum over k) 



_9 
dp 1 



p=0 



-^(^Tv{5(0) 75 A^(0) 7 fe ^(0) 75 }) 

x,y,z 



and we have considered the ratio 

0,0) 



R v (t) 



M 



(35) 



(36) 



p=0 



For large time separations Ry(t) is expected to be equal 
to the meson mass and, by setting the quark masses equal 
to the ones used in the previous section, it is possible to 
compare Ry(t) with mpp(t) and Rap^)- The compar- 
ison is shown in FIG. |2j As it can be seen the three 
definitions of effective meson mass agree within the sta- 
tistical errors in the middle of the lattice where the large 
time separation condition is satisfied also for Ry(t). 

We now come to the phenomenological relevant situa- 
tion in which the masses of the two mesons are different. 
In this case the determination of the two independent 
form factors can be achieved by solving the following lin- 
ear system of equations, 

<Pf|V°|Pj> = (Ej + E F )f + (q 2 ) + (Ej - E F )f_(q 2 ) , 



(PF\V k \P!) = (p k I+ PF)U(q 2 ) + (Pi-PF)f-(q 2 )- 

(37) 

The previous system becomes critical at the maximal 
value of <? 2 , i.e. in the zero recoil situation qj^ = 
(Mj — Mp) 2 . In this case the first equation can be used 
to obtain a single combination of /+ and /_ that is pro- 
portional to h+ defined here below (see ref. [3]) 



h±(q 2 ) 



(Mj + Mp) 

2y / M I M F 



f±(q 2 ) 



(Mj - Mp) 
lyJMjMp 



U(q 2 )- 



(38) 



The previous parametrization is widely used when the 
initial and final hadrons are both heavy-light pseu- 
doscalar mesons (for example in the B — >• Div case) and 



7 



the form factors are considered as functions of 



Pi -Pf 



MjM p 



(39) 



The zero recoil point at c^ M corresponds to uj = 1. 

By applying our method it is possible to obtain both 
the form factors directly at qj^. To this end it is useful 
to consider (by neglecting sub-leading exponentials) 



= I 



p=0 



\f^+f ](ah) Gl ° F e - M F{T/2-t)-M lt 
4MjM F 



g^ C iF(t^p) 



i>=o 



-\f-L- f \(a\A Gl ° F e -M F (T/2-t)-M lt 



(40) 



In previous expressions one of the quark lines has been 
drawn with a different color (green) in order to represent 
a quark of different mass with respect to the others (in 
particular in the numerical analysis we have set k green = 
0.13620 to be compared with k sea = 0.13590). By relying 
on the hermiticity of the vector current we also get 



c k FI (t,6,p) 



p=0 



- [fj.4- f MaiA Gl ° F e -Mi(T/2-t)-M F t 



Qpk 



c k FI (t,p,6) 



= I 



p=0 



\f+- f \(a\A Gl ° F e -M I (T/2-t)-M F t 



(41) 



By using these relations, together with the corresponding 
ones in the case of coinciding initial and final external 




FIG. 3: In the top panel we plot the effective form factor 
/+(<Zm) while in the bottom panel we show the effective form 
factor f-(q 2 M ) for which we get a good plateaux in the middle 
of the lattice where the large time separation condition is 
satisfied. 



states we get 

N n {t) = 
N FF (t) = 



AMj 



-M,T/2 



r 2 

-M F T/2 



(42) 



8 



The effective form factors are then obtained from 
R+(t)+R-(t) 



U(Qm) 



/-(<&) 



R+(t)-R-(t) 



(43) 



The combinations of correlators given in previous re- 
lations is particularly effective in minimizing statistical 
fluctuations that cancel between the different factors. 
In particular, when the initial and final meson masses 
are equal, the normalization conditions /+(0) = 1 and 
/-(0) = are identically verified for each gauge config- 
uration of the given ensemble. To give an idea of the 
plateaux that can be obtained by using this method we 
show in FIG. [3] the effective form factors extracted from 
our data. 



C. The vacuum polarization tensor 

In this section we show how the momentum expansion 
can profitably be used in order to calculate the hadronic 
vacuum polarization tensor, a quantity which is needed 
in order to predict the leading hadronic contribution to 
the anomalous magnetic moment of the muon. 

The hadronic vacuum polarization tensor is calculated 
on the lattice by considering the integrated two-point 
correlator of the electromagnetic currents of the quarks 



C^{p) 



1 



My . x+ , /2 . N2)(ytm{x)VM) 

2 ^^"in/Jl 



(TL 3 ) 2 

(S^f - p»p")Il(p 2 ) 



(44) 



where p^ = 2 sin(p /x /2). Given our choice of lattice Dirac 
operator, we have used the point-split vector currents 
in order to define the electromagnetic currents of the 
quarks, 



(45) 



We have chosen the momentum routing in which the ex- 
ternal momentum p flows through one of the two quark 
propagators. With this choice the point-split vector cur- 
rents at the vertices connect two fermion lines with differ- 
ent momentum and, consequently, bring the dependence 
upon p/2. 

The leading hadronic contribution to the g — 2 of the 
muon is obtained by extracting the scalar form factor 
II (p 2 ) from eq. (44) and by considering the difference 
H(p 2 ) — n(0). The subtracted form factor is introduced 
in order to cancel divergent contributions (see for exam- 
ple refs. [T2HT6] for dedicated works) and, for this reason, 
it is very important to have an accurate determination 
of 11(0). By using standard techniques, 11(0) can only be 
obtained by extrapolating the data obtained at p 2 > 0. 



These extrapolations unavoidably introduce systematic 
errors (for a recent discussion of this point see [8 ). In 
the following we shall show how 11(0) can be computed 
on the lattice directly and without the need of any ex- 
trapolation. 

As for the other quantities computed in this paper, 
the results for H(p 2 ) have been obtained with limited 
statistics, on a single volume, at fixed quark masses, etc. 
in order to demonstrate the effectiveness of the proposed 
procedure. For this reason we shall not use our results to 
make a prediction for the muon g — 2. This will be the 
subject of future work. 

After fermion integration and Wick contractions, 
the correlator C^ u (p) has both fermion-connected and 
fermion-disconnected contributions. If, for example, we 
limit ourself to the Nf = 2 case in the isospin symmetric 
theory we have 

C^(p) = (e„ + e d f <)0-( e '+ e l) O ■ W 

where e u = 2/3 and = — 1/3. In the following we 
shall concentrate on the connected part. Our method 
applies also to the disconnected part but we have not 
calculated the corresponding correlators and we shall not 
discuss these contributions in this work. Actually, many 
of the lattice collaborations involved in the computation 
of the hadronic vacuum polarization tensor [T2HT6] have 
neglected (or just attempted an estimate of) disconnected 
contributions. The choice, as in our case, is due to the 
big computational effort required for the calculation and 
not because these are expected (at least in general) to be 
negligible. 

In the following we call C^ v the connected contribu- 
tion corresponding to a single quark without multiplying 
it for the corresponding electric charge. In particular, 
we consider the correlators with different spatial indices 
(/i = k^v = h 7^ k) because these are proportional to 
the spatial momenta and because they do not require 
additional contributions in order to satisfy gauge Ward 
identities. Indeed, the lattice photon self energy with 
equal indices includes also the tadpole graph, where two 
photons couple to the fermion line in a single vertex (the 
current coupled to two photons is the one that we have 
previously indicated as "tadpole current"). We checked 
anyway, on limited statistics, that the gauge Ward iden- 
tities E^C^ = Y,vVvC^ v {p) = are satisfied. 

First, fixed \i = 1 and v = 2, we have computed the 
integrated correlation at p 1 > and p 2 > and divided 
it by the momenta, 



U(p 2 > 0) = 



c 12 ( P ) 

P 1 P 2 



1 



(tl 3 )' 



Tr {%, x; U\rb(x,m)S[x, V, U, X p }T 2 v (y,p/2)}) ■ 

(47) 

Then we have applied the rules discussed in the previous 
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i n(0) 
i n(p 2 >o) 



eq. pb) 



-0.1 0.1 0.2 0.3 0.4 0.5 0.6 

(*P) 2 

FIG. 4: Black points correspond to the calculation of 
n(p 2 ) p erformed by using standard techniques (according to 
eq. (47)). The red point correspond to 11(0) calculated di- 
rectly on the lattice (according to eq. (48 )). Data are in lattice 
units. 



sections to define the second mixed derivative, acting on 
propagators and vertices and evaluated at zero momen- 
tum, according to 



n(o) = 



d 2 C 12 {p) 
dp\dp2 



lp 2 =0 




(50) 



In FIG. [4] we show our results. The black points corre- 
spond to H(p 2 ) obtained from eq. (47) and, as expected, 
tend to be noisy for small values of p\ The red point cor- 
respond to 11(0) calculated directly on the lattice accord- 
ing to eq. ( 48]) . The uncertainty on 11(0) is, concerning 
our data obtained with limited statistics (~ 150 gauge 
configurations), comparable to the uncertainty that can 
be obtained at (ap) 2 ~ 0.05 but, coming from a direct 
calculation, it does not need to be corrected for system- 
atic errors due to extrapolations and, important to note, 
it scales with the statistics. 



n(o) 



d 2 C lz (p) 
dpidp 2 



(TL 3 ) 2 

p 2 =0 x,y 



VI. CONCLUSIONS 



Tr 



STl 



d 2 S 



-Tr 



v dpidp 2 V 

•2 



-Tr [ST^ST|] 



T d P[ 



-Tr 



dpi 



(48) 



where, for sake of brevity, we have dropped position ar- 
guments and we have used the relations 



dT k v (x,p/2) 
dpk 

dT k T {x,p/2) 
dpk 



T k T (x,p/2) , 
T k v (x,p/2) , 



(49) 



to obtain the derivative of the vertices (see section [IV] 
above). Note that in the previous expressions the factor 
1/2 appears because the currents here depend upon p/2 
and not upon p. 



To the lattice definition of 11(0), eq. (48), it can be 



given the following graphical representation (see also 



The method discussed in this letter allows the direct 
calculation on the lattice of the derivatives of correlators 
with respect to external momenta. We have described 
the method and checked its validity for several correlation 
functions. 

In particular, we have derived expressions to be used 
in order to compute both form factors parametrizing 
semileptonic decays of pseudoscalar mesons into other 
pseudoscalar mesons, directly at zero recoil. These rela- 
tions, checked numerically in this paper, may have many 
important phenomenological applications, for example in 
the calculation of B Dtv differential decay rate with- 
out excluding the t = r case, etc. Similar relations can 
be easily derived along the lines discussed in this paper 
for the pseudoscalar- to- vector channels. 

Another important application of the method concerns 
the direct calculation of the scalar form factor II (p 2 ) 
needed to predict the hadronic vacuum polarization con- 
tribution to the muon anomalous magnetic moment. In 
particular our method, as explicitly shown in the paper, 
allows to obtain 11(0) without the need of any extrapola- 
tion and related systematic uncertainties. 

The method is general, easy to implement, and we are 
pretty confident it will find several applications in the 
field of lattice phenomenology. 
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